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ABSTRACT 



Using the spectra of a sample of ~28,000 nearby obscured active galaxies 
from Data Release 7 of the Sloan Digital Sky Survey (SDSS), we probe the con- 
nection between AGN activity and star formation over a range of radial scales 
in the host galaxy. We use the extinction-corrected luminosity of the [OIII] 5007 
A line as a proxy of intrinsic AGN power and supermassive black hole (SMBH) 
accretion rate. The star formation rates (SFRs) are taken from the MPA-JHU 
value-added catalog and are measured through the 3" SDSS aperture. We con- 
struct matched samples of galaxies covering a range in redshifts. With increasing 
redshift, the projected aperture size encompasses increasing amounts of the host 
galaxy. This allows us to trace the radial distribution of star-formation as a 
function of AGN luminosity. We find that the star formation becomes more cen- 
trally concentrated with increasing AGN luminosity and Eddington ratio. This 
implies that such circumnuclear star formation is associated with AGN activ- 
ity, and that it increasingly dominates over omnipresent disk star formation at 
higher AGN luminosities, placing critical constraints on theoretical models that 
link host galaxy star formation and SMBH fueling. We parametrize this relation- 
ship and find that the star formation on radial scales <1.7 kpc, when including a 
constant disk component, has a sub-linear dependence on SMBH accretion rate: 
SFR oc M ' 36 , suggesting that angular momentum transfer through the disk 
limits accretion efficiency rather than the supply from stellar mass loss. 



1. Introduction 



Supermassive black holes (SMBHs) are intrinsically linked to their parent galaxies. 
As the black hole grows via accretion in the Active Galactic Nucleus (AGN) phase, the 
bulge of the host galaxy builds up from star formation, as evide nced by the correlation 



between SMBH mass and velocity dispersion of the host galaxy ( 


Kormendv & Richstone 


1995; 


Maeorrian et al. 


1998; 


Ferrarese & Merritt 


2000; 


Gebhardt et al. 


2000; 


Tremaine et al. 
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20021 ; lHaring fc Rixl 12004 ) . Although the evolution of SMBH accretion and star forma- 
tion are surely intertwined, the mechanisms that couple these processes are not defini- 
tively known. Gas can be driven to the potential well of the black hole by galaxy mergers, 



which compress gas and thereby trigger host galaxy star formation (e.g., [Sanders et al 



Kauffmann Sz Haehnelt 2000 ; Di Matteo et al. 2005 ; Springel et al. 2005 ; [Hopkins et al. I l2005h 



1988; 



Alternatively , secular mechanisms such as torques from galactic bars, bars-within-bars 



stochastic accretion of cold gas (IHopkins fc Hernquistl 



f lNorman fc Scovilld 11988 : ICiotti fc Ostriker 



spiral arms (IShlosman et al.l Il989t iGarcia-Burillo et al.l 120051 ; IHopkins fc Quataertl |2010| ) 



1997 



2007 



20061) or mass loss from stellar winds 



Kauffmann fc Heckmanl |2009| ) can 



contribute significantly to both processes. Feedback is another key ingredient in the star for- 
mation and AGN connection, where AGN winds can blow out circumnuclear gas, suppressing 



star f ormation and regulat i ng SMBH and bulge g r owth flSilk fc Rees 



1998 



Granato et al. 2004; Di Matteo et al. 2005 



Wvithe fc Loeb 



19981 
2003 



Haehnelt et al. 



Somerville et al. 



20081 ). Observational studies of local optically obscured AGN (Type 2 Seyferts, Sy2s) have 
indicated that star formation processes can both feed the black hole, via s tellar mass loss, 
and limit accretion efficiency through supernova feedback (jWild et al.ll2010l ). 



Seyfert 2 galaxies are prime candidates for studying the observational link between 
SMBH growth and host galaxy star formation. As the accretion disk is enshrouded by 
circumnuclear dust obscuration, optical signatures of star formation in the host galaxy as well 
as AGN fueling can be investigated. Previous studies of Sy2s have revealed that a significant 
fraction live in host galaxi es that are undergoing (or have recently experienced) significan t 
amounts of star formation fjKauffmann et aL 2003a ; Cid Fernandes et al. 2004 ; Netzer 2009 ). 
Usin g a sample of over 20,0 AGN from the S l oan Di gital Sky Survey (SDSS) Data Release 
One fIStoughton et alJl2002h , [ Kauffmann et al.l (j2003al ) demonstrated that the trend between 
star formation activity and AGN accretion wa s not confined to the n uclear regions of the 
host galaxies for luminous AGN. Subsequently, iKauffmann et al.l (120071 ) used a combination 
of SDSS spectra and GALEX UV images to show that the link between star- formation and 
black hole growth was much more direct for star formation in the ce ntral few kpc than for star 
forma tion in the global galactic disk. In agreement with this result, iDiamond-Stanic fc Rieke 
(120 121) studied 74 galaxies from the revised Sh ap ley- Ames Seyfert sample f lShapley fc Ames 



1932 



Maiolino fc Riekd Il995t iHo et al.l Il997f ) and found that this connection is strongest 



on sub-kpc scales, with a much weaker correlation when considering larger scale (i.e., >1 
kpc) star formation. These place important constraints on the mechanisms by which star 
formation and black hole growth are linked. 

However, these studies suffer from non-uniform methodology for measuring star forma- 
tion rates and/or small sample sizes. To address these limitations, in this work we test the 
relationship between star formation and SMBH accretion on various galactic scales using a 
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sample of ~28,000 type 2 AGN and composite galaxies culled from SDSS. All quantities are 
measured in a uniform way through the 3" spectroscopic fiber. These quantities parametrize 
the AGN luminosity, SMBH accretion rate, and star formation rate (SFR). We test whether 
the connection between AGN fueling and starburst activity varies as the projected size of 
the SDSS aperture sequentially covers greater parts of the host galaxy, and whether this 
relationship is a function of intrinsic AGN luminosity. 



2. Data Analysis 



Our data were drawn from the Main Galaxy Sample of SDSS Data Release 7 (DR7 
Strauss et al.l |2002| ; lAbazajian et al.l 120091 ) where we utilized the MPA-JHU value added 
catalogs^ We isolated s ources categorized a s "AGN" and "Composite" according to th e 
diagnostic BPT diagram ( iBaldwin et al.l 1 198 It iKewley et al.ll2001t iKauffmann et al.ll2003al ); 
we note these classifications were made automatically by the MPA-JHU pipeline. We ex- 
clude galaxies classified as LINERs. We imposed several quality control cuts: the Ha, H/3, 
[OIII] 5007A and [Nil] 6584 A emission lines had to be significant above the 5<r level and the 
Ha equivalent width had to exceed 3 A removing galaxies where the ionizing radiation may 
be due to low- mass evolved stars and not nuclear accretion ( ICid Fernandes et al.ll201ll ). We 
also restrict the sample to the Petrosian magnitude range 15.5 < m r < 17.8 (the limiting 
magnitude of the Main Galaxy Sample) and remove outliers in BPT parameter space, focus- 
ing between -1.0 < log([NII]/Ha) < 1.0 and -1.0 < .og([OIII1/H/.) < 2 . Th e objects in our 



sampl e span 0.01 < z < 0.32. As shown by IKauffmann et al.l (]2003al ) and iReichard et al. 



(120091 ). these AGN host galaxies are typically structurally normal galaxies of intermediate 
Hubble type (average SDSS concentration index of 2.61±0.33). 



2.1. AGN Activity 



We use the luminosity of the [OIII] line {h\oni\) as a proxy for intrinsic AGN luminosity 
since it is primarily ionized by accretion disk photons. Since [OIII] forms in the AGN narrow 
line region, it is not affected strongly by dust obscuration from the circum nuclear torus, pro 



viding a good estimate of the intrinsic AGN lumin osity in SDSS galaxies (IKauffmann et al. 



2003a; Heckman et al. 



2001; 



LaMassa et al.ll2010l ). For composite galaxies, i.e., those sys- 



tems between the Kauffmann et al. (l2003alj empirical boundary between star-forming and 
active galaxies and the IKewley et al. ( 2001 ) theoretical maximum starburst line on the BPT 



1 http://www.sdss3.org/dr9/algorithms/galaxy_mpa_jhu.php 
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diagra m, we estimated t he AGN contribution to the [OIII] emission following the prescrip- 
tion of I Wild et al.l (120 10[ ). This method estimates the star formation contributions to L\oun 
based on the location of the source in the BPT diagram. We only keep composite galaxies if 
the AGN component of the [OIII] line is at least 30% of the total [OIII] emission, removing 
very weak AGN. In total, we have 27952 galaxies, of which 15508 are composites and 12444 
are Sy2s; they are shown on the BPT diagram in Figured! Finally, we correct ~L[oni] for dust 
obscuration using the observed Balmer decrement (Ha/H/3), assuming an in trinsic Ha/H/3 



value of 3.1, and the standard R(V) = A(V)/E(B — V)=3.1 extinction curve (ICardelli et al. 



19891 ). obtaining L[oni], 



The SMB H mas (Mbh) goes a s the host galaxy stellar velocity dispersion (cr) to the 
fourth power (jTremaine et al. 2002 ). To parametrize the SMBH accretion rate, we there- 



fore define a pseudo-Eddington parameter, L 



/a 4 , for the galaxies where a velocity 



[OIII] ,corr 

dispersion was measured. In the following, we search for trends between AGN luminosity 
(L[oni],corr) and SFR, and between the Eddington parameter (L[oni],corr/o' 4: ) and specific 
SFR (sSFR, the SFR normalized by galaxy mass). 



2.2. Star Formation 



We focus on the SFR and sSFR measu red through the 3" SPSS spectroscopic fiber 
(SFR/jter, sSFR/jb er ), calculated as outlined in lBrinchmann et al.l (120041 ) . As nuclear activity 
contributes to Ha luminosity, which is used to calibrate the SFR in quiescent galaxies, these 
SFRs are instead c alculated using the break in the optical spectrum at 4000 A (D n (4000), 
Balogh et al.lll999l ). which is a measure of the age of the host galaxy stellar population: a 
shallower break is associated with a younger stellar population. Using t he relationship be 
tween Dn(4000) and the Ha calibrated sSFR/j{, er for star-forming galaxies. iBrinchmann et al. 



(120041 ) constructed probability distribution functions from which sSFR/j6 er . were estimated 
for active galaxies; SFRyj ber is then calculated by multiplying sSFRjj 6er , by the galaxy mass 
within the fiber. In this analysis, we use the value-added catalog provided median estimates 
of the SFR/jf, er and sSFR/jf, er probability distribution functions. 



2.3. z/z max as a Distance Proxy 

Our goal is to determine the relationship between black hole growth and star formation 
over a range of physical scales inside the galaxy. To ensure that we are comparing host 
galaxies with similar properties at all distances, and isolating the physical size of the pro- 
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j ected SPSS fiber apertu re as the varying parameter, we introduce the z/z max parameter 



( iKauffmann et al.ll2003al Jbl). where z max is the redshift at which the galaxy falls out of the 
apparent magnitude limit of the sample (m r = 17.8). Absent significant cosmic evolution 
(see below), the properties of a population of galaxies selected from SDSS will be invariant 
with z/z m ax- Normalizing the redshift of each galaxy by z max corrects for the effect that 
more luminous and massive galaxies are visible out to greater distances. To calculate z ma x, 
we first calculate the absolute r magnitude (M r ) using: 

M R = m r - DM(z) - k(z), (1) 

where DM(z) is the distance modulus and k(z) is the k-correction. We adopted a 
cosmology of H = 70 km s" 1 Mpc" 1 , Q M = 0.27 and A=0.73. 



We determined k-corrections for each galaxy using kcorrect version 4.2 ( IBlanton fc Roweis 



20071 ). Using we then found the redshift, z max , at which each galaxy would have an 
apparent magnitude of 17.8. As indicated in the legends of Figure [2J z/z max is proportional 
to redshift. 



2.4. The Effects of Cosmic Evolution Sc Sample Completeness 

With our approach of using galaxies over a range of redshifts to study the radial distri- 
bution of star-formation in AGN, there are two issues that must be addressed. Firstly, since 
we are confining our study to the relatively local universe, we do not expect the number 
density of AGN to evolve with redshift. Our sample binned in z/z max only span a range in 
median redshift from 0.06 to 0.13. Over this range, the increase in the rate of star-formation 
and black hole growth per co-moving volume element (e.g. Marconi et al. 2004) is only a 
factor of ~ 1.2 (0.08 dex). Secondly, low-luminosity AGN will be increasingly difficult to 
recognize in more distant galaxies since the diagnostic emission-lines will become fainter and 
more diluted by emission due to host galaxy star formation. 

To address these two issues, we performed a V/V max test for our sample, where V is 
the cosmic volume probed by an individual galaxy sample and V max is the maximum volume 
out to which that galaxy would be observable before dropping below the sensitivity limit of 
the survey. Values for the median of a V/V max distribution of significantly less than (greater 
than) 0.5 signals incompleteness (cosmic evolution). In TableHJ we list the median V/V max for 
various AGN luminosity and Eddington parameter bins. In the following, we work within the 
range where V/V ma x = 0.5 i 0.1, i.e., L[oiii],corr 

> 10 40 erg s" 1 and L [0 iii],corr/^ > 3 x 10 31 
erg s 3 km -4 . We also find that the median estimate of the total stellar mass and the velocity 
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dispersion does not evolve with z/z max in these AGN luminosity and Eddington parameter 
bins. Thus, we are comparing similar galaxies in aggregate over increasing z/z max values. 

Adopting the bolometric c orrection for the extinction-corrected [OIII] luminosity from 
Kauffmann fc Heckmanl (120091 ) and the conversion from a to Mbh from iTremaine et al. 
( 120021 ) . our sample covers ranges from roughly 10 9 ' 2 to 10 12 ' 2 L Q in bolometric luminosity 
and from roughly 10~ 2 ' 8 to 1 in L/LiEdd- Over these broad ranges we are reasonably complete 
and the effects of cosmic evolution are modest. 



3. Results 

Separating galaxies into bins of z/z max , we plot SFR/jft er as a function of AGN luminosity 
(L[onr\,ccrrr) an d sSFRyj feer as a function of the Eddington parameter (L[ ///] corr /cr 4 ), and 
the associated dispersions, in Figures [2J Within each z/z max bin, we further bin in AGN 
luminosity and Eddington parameter, including at least 20 sources in each bin. We then 
calculate the median SFR and sSFR in these bins. The dispersion (disp) is the associated 
la error for a Gaussian distribution, i.e., half the difference between the 83rd and 17th 
percentiles of the SFR/jb er and sSFR/j& er distributions in each bin. 

As z/z max increases, the 3" spectroscopic fiber covers larger physical scales of the host 
galaxy. The median galaxy redshift ranges from ~0.06 to ~0.13 between the lowest and 
highest z/z max bin, corresponding to physical aperture radii of 1.7 - 3.5 kpc. As can be seen 
in Figure [21 the SFR (sSFR) increases slowly but systematically with the AGN luminosity 
(Eddington ratio), within a given z/z max bin. It is also clear that there are systematic 
differences between these relations as a function of SDSS aperture size (z/z max ). At low to 
moderate AGN luminosities there is a trend for the SFR to grow for a given AGN luminosity 
as the projected aperture size increases, although this trend is not clear at higher AGN 
luminosities. Moreover, the slope of the relationship between sSFR and the Eddington-ratio 
proxy flattens as z/z max increases. 

To see these trends with distance more clearly, in the left panels of Figure [3] we plot the 
star formation rate (top panel) and sSFR (bottom panel) as functions of z/z max , and bin the 
data by the AGN luminosity (top) and Eddington ratio (bottom). These plots demonstrate 
that the radial variations of SFR and sSFR are dependent on the AGN luminosity and 
Eddington ratio, respectively. To visualize and quantify how the radial distribution of SFR 
and sSFR vary as a function of increasing AGN luminosity and Eddington ratio, we have 
calculated the slope, a, of the correlation between 



(s)SFR oc (z/z max ) a , 



(2) 
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in each AGN luminosity and Eddington parameter bin. For the linear regression fit, we 
weighted the data points by the error in each bin, taken as disp/y/N, where disp is the 
dispersion in the distribution within the bin (half the difference between the 17% and 83% 
percentiles) and N is the number of sources per bin. We note that the lowest z/z max bins 
have the smallest number of sources, which results in the associated errors being larger. 

The results in the right panels of Figure [3] show the same systematic trend. As the AGN 
luminosity increases, the increase in star-formation inside the SDSS aperture becomes smaller 
with greater projected aperture size. At lower AGN luminosities, the relationship between 
SFR and increasing projected size of the host galaxy is roughly linear, showing that star 
formation is distributed over the whole range in radial scales we probe in these host galaxies. 
However, at higher AGN luminosities, star formation is predominantly circumnuclear: we do 
not see evidence that the SFR increases significantly as the size of region we observe increases 
from 1.7 to 3.5 kpc (the physical size of the fiber radius at the median redshifts from lowest 
to highest z/z max bins). Similarly, at high Eddington rates, the negative slope between 
sSFR and projected aperture size further suggests dominant circumnuclear star formation: 
as the size of the observed region grows, stellar mass is added, but is not compensated by 
an increase in star formation. 



4. Discussion 

These results indicate that circumnuclear star formation associated with AGN activity 
(SFK nuc i eus ) only dominates over omnipresent disk star formation (SFR^sfc) when the AGN 
is very luminous. During such a stage, host galaxy star formation will appear to be centrally 
concentrated. Let us assume 

SFR 

fiber — S F ' Rnucleus ~\~ S F Rdiski (3) 

The star formation in the nucleus is associated with AGN fueling, parametrized by 

SFR nudeus = a(M Q /yr) x {L [OIII]jCorr /10 42 erg/ S y . (4) 

We use SFR/jb er from the lowest AGN luminosity bin, 40.0 < log {J->[oiii],corr) < 40.5 dex, 
as the omnipresent SFR^fe. We solve equation [3] for all log (L[oiii},corr) > 40.5 dex bins 
in Figure [3] (a) simultaneously by minimizing x 2 . Here we take the error on each bin to 
be disp/y/N, where disp is half the difference between the 83rd and 17th percentile and N 
refers to the number of sources in the bin. We find a = 0.44 ± 0.02 and = 0.36 ± 0.04, 
and plot these fits in Figure [3] (a); the errors represent the 95% confidence interval, where 
A% 2 =3.84. These results indicate that star formation associated with SMBH fueling rises 
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sub-linearly with AGN luminosity. The fit could be improved with a more complex model 
in which both the disk and circumnuclear star formation rates are allowed to have separate 
functional dependencies on the AGN luminosity, but such complexity is not really warranted. 
Our simple toy model is only meant to show that the data are consistent with a stronger 
link between the AGN luminosity and the amount of central star-formation than with larger 
scale star-formation. 

What can be the physical link between these two processes? Circumnuclear star- 
formation in luminous AGN can originate from major mergers, as simulations have demon- 
strated that torques produced by galaxy-galaxy interactions fuel SMBHs while triggering star 



form ation in t he galactic cen t er (|Di Matteo et all 120051 ; ISpringel et al.l 120051 ; iHopkins et al. 



20051 ). Indeed, iTreister et al.l (120 121 ) have demonstrated that mergers are responsible for trig- 
gering the most luminous AGN (Ly > 10 45 erg s _1 ). However, the majority of the highest 
luminosity AGN in our sample falls below this bolometric luminosity threshold. Additionally, 
the lopsidedness of SDSS Sy2 host galaxies within the Eddington range we probe (~ 10~ 3 
to one) does not drastically increase, suggesting that mergers are n ot the primary mecha - 
nism for coup ling AGN and star form ation activity in this sample ( iReichard et al.ll2009l ). 
Alternatively, ICiotti fc Ostrikerl ( 120071 ) have shown that recycled gas from dying stars near 
the galactic centers of elliptical galaxies becomes radiatively unstable, causing circumnuclear 
star formation while also feeding the SMBH. As instabilities caused by bars, bars-in-bars, 
spiral arms, etc. produce galactic scale star formation while funneling gas to the SMBH 
(IHopkins fc Quataertll2010l ). stellar mass loss seems a more likely mechanism for coupling 
AGN activity and circumnuclear star formation in moderate luminosity Sy2 galaxies. 



Conclusions 



We find that the AGN-star formation connecti on is largely driven by c ircu mnuclear star- 

burst activity in local Type 2 AGN, consistent with lKauffmann et al.l (120071 ) and lDiamond-Stanic fc Rieke 
(l2012h . 



The size-scale of star formation in local active galaxies appears to depend on AGN 
luminosity: it becomes increasing compact in more luminous AGN. More specifically, the 
lack of an increase in the amount of SFR as the size of the region probed by the SDSS 
spectra increases implies that star formation is concentrated in radii below <1.7 kpc in 
these luminous systems. This constrains the mechanism that couples AGN activity and host 
galaxy star formation, favoring theories that link these two processes at circumnuclear rather 
than galactic scales (i.e., stellar mass loss rather than instabilities caused from bars, etc.). 
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Our results are consistent with the hypothesis that circumnuclear star formation is 
associated with AGN activity and thus dominates over omnipresent galactic disk star for- 
mation when the AGN becomes luminous. The dependence of the SFR on SMBH fueling is 
sub-linear: on radial scales <1.7 kpc, SFR oc M 36 . This sub-linearity may indicate that 
accretion is a "demand" rather than "supply" driven process. Angular momentum trans- 
fer through the disk may be relatively inefficient, limiting the accretion rate, regardless of 
the amount of fuel provided via stellar mass loss. Remaining material could potentially be 
ejected in winds and/or jets, or perhaps accumulate in the obscuring "torus." 
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Fig. 1. — BPT diagram for galaxies used in our analysis. The lower dashed l ine represents the 
empiric al boundary between star forming galaxies and active systems from iKauffmann et al. 
(I2003af) while the up per dotted line delineates the maximum theoretical starburst line from 
Kewley et al.l (120011 ) . There are 15508 composite galaxies (red data points, between the two 
curves) and 12444 Sy2 galaxies (blue data points). These composite systems represent the 
subset of composites in which the AGN contribution to [OIII] emission is estimated at > 
30% of the total. 
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Table 1. Median V/V max for AGN Luminosity and Eddington Parameter Bins 



Log (L[ /77] jCorr ) V/V max Log (L[ ///] iCorr /<7 4 ) V/Vmax 

erg s _1 erg s 3 km~ 4 



38.5 - 39.0 


0.25 


30.0 - 30.5 


0.36 


39.0 - 39.5 


0.33 


30.5 - 31.0 


0.32 


39.5 - 40.0 


0.38 


31.0 - 31.5 


0.37 


40.0 - 40.5 


0.44 


31.5 - 32.0 


0.43 


40.5 - 41.0 


0.49 
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0.46 


41.0 - 41.5 


0.51 


32.5 - 33.0 


0.50 
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0.55 
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0.53 
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0.59 
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Fig. 2. — (a) SFR/jber as a function of AGN luminosity and (c) sSFRyjb er as a function of Ed- 
dington parameter in increasing z/z max bins, (b,d) associated dispersions and (c,e) legends. 
The legends give the median redshift of the galaxies in each bin: with increasing redshift, 
the projected physical aperture size grows, encompassing greater amounts of the host galaxy. 
The relationship between AGN luminosity and star formation shows a dependence on the 
amount of the host galaxy that is sampled at lower AGN luminosities, but not at higher 
ones. The slope of the relation between sSFR and Eddington parameter becomes flatter 
as the aperture size grows, meaning the star formation in more luminous AGN is centrally 
concentrated. 
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Fig. 3. — (a) Log (SFRy ifeer ) as a function of z/z max in increasing L[oni],carr bins and (b) 
the best fit slope in each bin. (c) Log (sSFR/j& er ) as a function of z/z max in increasing 
L[oi"//],corr/c" 4 bins and (d) the best fit slope in each bin. At higher AGN luminosities and 
accretion rates, galactic star formation becomes more centrally concentrated. The lines 
in (a) correspond to the best- fit parameters from SFRj iber = SFR nuc ; ear + SFR disfc , where 
SFR nudear = a (M /yr) x(L[ O /ii] )CorT /10 42 erg/s)' 3 and SFR disfe is SFR /ifeer associated with 
the lowest AGN luminosity bin (i.e., 40.0 < log(L[ m] iCor r < 40.5 dex). From a reduced x 2 
minimization, we find a=0.44±0.02 and /3=0.36±0.04. 
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